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We perform first-principles calculations of electronic structure and optical properties for UO2 
and PuC>2 based on the density functional theory using the generalized gradient approximation 
(GGA)+{7 scheme. The main features in orbital-resolved partial density of states for occupied / 
and p orbitals, unoccupied d orbitals, and related gaps are well reproduced compared to experimental 
observations. Based on the satisfactory ground-state electronic structure calculations, the dynamical 
dielectric function and related optical spectra, i.e., the reflectivity, adsorption coefficient, energy-loss, 
and refractive index spectrum, are obtained. These results are consistent well with the attainable 
experiments. 



PACS numbers: 78.20.-e, 77.22. Ch, 71.20.-b 
I. INTRODUCTION 

Actinide dioxides (An02) have been attracted lots of 
attention due to their rich physical phenomena charac- 
terized by the complex nature of 5/ electrons. Many 
experimental and theoretical works have been devoted to 
investigating the thermodynamical, electronic structural, 
and defect properties of AnC>2 systems. Taking UO2 and 
PuC>2 for example, their insulating ground states have 
been established experimentally Ju 0] and successfully 
predicted theoretically 0, 0, S W- When referring to 
insulators or semiconductors, one basic physical quan- 
tity of interest is their band gaps. If the band gap of 
UO2 or PuC>2 can be comparable to semiconductors, one 
idea may occur to us that whether they can be applied 
extensively in the electronic and optoelectronic devices 
like semiconductors (Si, GaAs, and ZnO) or not. Re- 
cently, Meek et al. discussed the electronic properties of 
uranium dioxide and revealed the potential performance 
advantages of uranium dioxide as compared to conven- 
tional semiconductor materials Q • Especially, the higher 
dielectric constant of UO2 makes it more suitable for 
making integrated circuits [7( . This may stimulate many 
studies of the optical properties for actinide dioxides in 
future. 

Optical adsorption and reflectance spectra of semicon- 
ductors have been studied for several decades both exper- 
imentally and theoretically, whereas, similar works per- 
formed on actinide dioxides is still very scarce although 
they are necessary not only from the viewpoint of ba- 
sic science but also from their technological importance 
in industries. Experimentally, Schoenes studied the in- 
cidence reflectivity of UO2 single crystals in the photon 
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energy range of 0.03-13 eV, from which the complex di- 
electric function e{oj) — £i(u>) + ie-i (lu) has been derived 
[8[. For PuC>2, to our knowledge, no experimental opti- 
cal data are available in literature. As for the theoretical 
investigations of optical spectrum of actinide dioxides, it 
is a great challenge to standard density functional the- 
ory that an accurate description of electronic structure 
for actinide oxides is hard to be achieved, which is indis- 
pensable to getting the correct optical spectrum. Con- 
ventional density functional schemes that apply the local 
density approximation (LDA) or the GGA underestimate 
the strong on-site Coulomb repulsion of the 5/ electrons 
and consequently fail to capture the correlation-driven 
localization. Therefore, the 5/ electrons in actinide ox- 
ides require special attention. One promising way to 
improve contemporary LDA and GGA approaches is to 
modify the intra-atomic Coulomb interaction through the 
so-called LDA+U or GGA+f/ approach, in which the 
underestimation of the intraband Coulomb interaction 
is corrected by the Hubbard U parameter |T3|. Re- 
cently, the electronic structures of UO2 and Pu02 are 
correctly reproduced using LDA+ U or GGA+ U calcu- 
lations [1, |j, _3, __[. Therefore, based on the good per- 
formance of LDA/GGA+ U approaches in describing the 
electronic structure of the systems containing 5/ elec- 
trons, it is encouraging to investigate the optical spectra 
of them. 

In this work, we used the GGA+ U scheme to study 
the static and frequency-dependent dynamical dielectric 
response functions for UO2 and Pu02- Our present cal- 
culated band gap E g and high-frequency dielectric con- 
stant £oo for UO2 are 2.3 eV and 5.53, which are in good 
agrement with the experimental values of about 2.1 eV 
and 5.1 observed in the optical spectra @, respectively. 
Furthermore, our calculated dielectric function e(u) ex- 
hibits the overall agreement with experimental result and 
the main peaks are well reproduced. The dielectric func- 
tion and the consequent optical spectra for Pu02 are also 
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calculated in the paper. In particular, the value of Eqo for 
PuC>2 is predicted to be 6.21, a little larger than that for 
UO2. Considering the satisfactory calculations for UO2, 
we expect our predicted optical behavior for PuC"2 can 
provide a useful reference for future experimental mea- 
surement. 



II. DETAILS OF CALCULATION 



following equation [15| : 
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where the indices c and v refer to conduction and valence 
band states respectively, and u c k is the cell periodic part 
of the wavefunctions at the fc-point k. The real part of 
the dielectric tensor is obtained by the usual Kramers- 
Kronig transformation 



Our electronic structural and optical calculations are 
performed using the projector-augmented wave (PAW) 
method of Blochl as implemented in the ab ini- 

tio total-energy and molecular-dynamics prog ram VASP 
(Vienna ab initio simulation program) [12j. PAW is 
an all-electron method that combines the accuracy of 
augmented-plane-wave methods with the efficiency of 
the pseudopotential approach. The PAW method is 
implemented in VASP with the frozen-core approxima- 
tion. The exchange-correlation functional is used GGA 
of Perdew-Burke-Ernzerhof (PBE) formalism [l3| . The 
5/, 6s, 6p, 6d and 7s electrons of U and Pu as well as 
the oxygen 2s and 2p electrons are explicitly treated 
as valence electrons. The electron wave function is ex- 
panded in plane waves up to a cutoff energy of 500 eV. 
For the Brillouin zone integration, the T centered 6x6x6 
grid is adopted. 144 bands are used to get the dynami- 
cal dielectric function e(uj) and a good convergence can 
be achieved. In order to perform the antiferromagnetic 
(AFM) phase calculations, we used the unit cell contain- 
ing 12 atoms. The strong on-site Coulomb repulsion 
among the localized 5/ electrons is described by using 
the formalism formulated by Dudarev et al. [14| . In 
this scheme, only the difference between the spherically 
averaged screened Coulomb energy U and the exchange 
energy J is important for the total LDA (GGA) energy 
functional. Thus, in the following we label them as one 
single effective parameter U for brevity. In our calcula- 
tion, we use J=0.51 and 0.75 eV for the exchange energies 
of U and Pu, respectively, and the effective Hubbard U 
are 4.0 and 3 eV, which are close to the values used in 
other previous work 0, Q . 

As for the optical spectra calculations, we adopt two 
different methods to determine the macroscopic static 
dielectric constants using different approximations [l5l |. 
One method is using a summation over conduction band 
states and the other is using the linear response the- 
ory (density functional theory). For the latter, only the 
static ion-clamped dielectric matrix can be obtained and 
a summation over empty conduction band states is not 
required, whereas the former can calculate the frequency- 
dependent dynamic dielectric function after the elec- 
tronic ground state has been obtained. The frequency- 
dependent imaginary part of the dielectric function is de- 
termined by a summation over empty states using the 
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where P denotes the principle value. 

The main optical spectra, such as the reflectivity R(lj), 
adsorption coefficient /(w), energy-loss spectrum L(u>), 
and refractive index n(uj), all can be obtained from the 
dynamical dielectric response functions s(lv). The ex- 
plicit expressions are given by 
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respectively. 

III. RESULTS AND DISCUSSIONS 



electronic structure and optical properties of 
U0 2 



Since the optical spectra are directly calculated from 
interband transitions, an accurate description of the elec- 
tronic structure is indispensable. The calculated orbital- 
resolved partial density of states (PDOS) for U 5/, U 6d 
and O 2p are shown in Fig. 1(a). The Fermi level is set 
to be zero. It is clearly shown that the valence bands are 
mainly contributed by U 5/ and O 2p orbitals. The peak 
near the Fermi level is mainly U 5/ with a little O 2p 
contribution, which bas been confirmed by the resonant 
photoemission[l6j]. The U 5/ valence band covers from 
to —1.6 eV, which is also consistent with the experi- 
mental observation that the occupied 5/ states in UO2 
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FIG. 1: (a) The projected orbital-resolved partial DOS for U 
6d, U 5/, and O 2p orbitals in antiferromagnetic UO2. The 
Fermi level is set to zero, (b) The dynamical dielectric func- 
tion e(w) = Ei(uj) + C2(w) as a function of the photon energy 
iv for UO2. The black and green lines represent our calculated 
real and imaginary parts of dielectric function e(uj), respec- 
tively, while the red dotted-line is experimental £2(0;). 
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FIG. 2: Calculated optical spectra for UO2, (a) the reflectivity 
R(io), (b) adsorption coefficient I(u>), (c) energy-loss L(lu), 
and (d) refractive index n(w). 



TABLE I: Ion clamped static macroscopic dielectric constants 
£00 of UO2 and PUO2 calculated using the PAW method and 
various approximations with various fc-points sampling: V in- 
dicates a grid centered at F point, whereas Monkhorst-Pack 
(MP) grids do not contain the F point. stands the number 
of irreducible fc-points of the Brillouin zone (IBZ) at specific 
fc-points sampling. e m i c indicates values neglecting local field 
effects, srpa includes local fields effects in the Hartree ap- 
proximation, and £dft includes local fields effects on the DFT 
level. e cond are values obtained by summation over conduc- 
tion band states, whereas £ LR are values obtained using linear 
response theory (density functional perturbation theory). 
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are located around 1.5 eV below the Fermi level with a 
band width of about 2.0 eV The O 2p valence band 
width is 4.0 eV from about —1.8 to —5.8 eV, in quali- 
tative agreement with the photoemission value of 5.0 eV 
from -3.0 to from -8.0 eV 

As for the unoccupied U 5/ and 6d orbitals, their accu- 
rate descriptions are also indispensable to the interband 
transitions, since electrons are excited from the occupied 
valence bands to the unoccupied bands during optical ex- 
citations. The 5/ and 6c? bands begin at about 2.3 and 
4 eV, respectively, which are consistent well with the re- 
sults of 2.6 and 5 eV obtained by hybrid DFT method 
Note that our calculated p — > d gap is 5.8 eV, which 
accords well with the Bremsstrahlung Isochromat Spec- 
troscopy (BIS) value of 5.0±0.4 eV [lj]. Overall, our 
calculated DOS agrees well the experimental spectra and 
other theoretical results. This supplies the safeguard for 
our following optical spectrum calculations. 

Due to the cubic symmetry of UO2, the dielec- 
tric tensor only has one independent component and 
£xx=£yy=£zz- Our calculated macroscopic dielectric con- 
stants £oo using different methods and approximations 
are collected in Table I. We find that well converged re- 
sults can be obtained by using the T-centered 6x6x6 
grid. Note that the value of £^ T should be compared to 
experiment. For UO2, the calculated Soc is 5.53, which 
agrees well with the experimental value of 5.1 @. 

As for the dynamical dielectric function, our calculated 
imaginary part £2(1^) and real part £\{lo) of the complex 
dielectric function e(oS) together with the corresponding 
experimental £2(1^) are showed in Fig. 1(b). The green 
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and black lines represent our calculated imaginary and 
real parts of the complex dielectric function b(lu), respec- 
tively, while the red dotted-line gives the experimental 
measurement @ of £2(w)- Our theoretical photon energy 
covers from to 30 eV, while the experimental [|| value 
covers from to 13 eV. According to our calculated DOS 
showed in Fig. 1(a), we suggest that in £2(w) the peaks 
(at 2.8 eV) below 3 eV should be assigned to the in- 
tra 5/ transitions. Notice that the unoccupied 6c? bands 
begin about at 4 eV, therefore, the 5/ — > 6c! transition 
energies should be larger than 4 eV. Kudin et al. also 
suggested that the stronger adsorption observed experi- 
mentally at ~5-6 eV could be assigned to the optically 
allowed 5/ — > 6d transitions [l7j . According to our cal- 
culated £2(w), four main peaks lie at about 5.0, 7.1, 9.8, 
and 11.8 eV, respectively. The shape of the calculated 
curve exhibits the same main features demonstrated by 
the experimental results [8j. Combined with the orbital- 
resolved PDOS shown in Fig. 1(a), we attribute the first 
two peaks in £2(1^) to be 5/ — * 6c! transitions, while the 
last two to be 2p — > 6d transitions. This is consistent 
well with the experimental assignment [l!| by Naegele 
et al., who attributed the peak around 3 eV in e-ziiS) to 
intra 5/ 2 transitions, while the peak structures above 5 
and 10 eV were ascribed to the / — > d and p — > d transi- 
tions, respectively. Another assignment was suggested by 
Schoenes according to their dielectric function deduced 
from the reflectivity measurement; they argued that the 
peaks near 3 and 6 eV correspond to / — > d transitions, 
and that the peaks near 8 and 11 eV are due to p — > 
d transitions [2(|. Herein, the assignment of / — ► d 
transition at 3 eV in Ref. [1, [2(| is not supported by our 
calculation. The cause is that in assigning the peak in 
Ei{u) at 3 eV, the energy distance between U occupied 
5f 2 and O 2p valence bands was overestimated in Ref. 
■8j, 20] to be as large as 4 eV, which is much larger than 
that directly determined by the photoemission measure- 
ments [16|, l2l|, |22[. On the contrary, according to our 
band-structure calculation, the occupied 5/ orbitals are 
locate at about 1.5 eV below the Fermi level and the O 2p 
bands widely covers from about -1.8 to -5.8 eV, which in- 
stead accords well with the experimental photoemission 
data [H, [HI, HH in UO2. Thus, as mentioned above, we 
suggest the structure in £2^) below 3 eV is caused by 
the intra 5/ transitions. 

Using expressions (3)-(7), the reflectivity R(uj), ad- 
sorption coefficient I (to), energy-loss L(u>) and refractive 
index n(u>) spectra are showed in Fig. 2. For reflectiv- 
ity R(u>) spectrum, there are four peaks locating at 4.8, 
7.6, 10.3, and 12.8 eV. The adsorption coefficient 
spectrum has the same trends. The origin of these peaks 
can also be explained as the peaks of the imaginary part 
£2(k>). Note that three similar peaks at 5.5, 8, 11.7 eV 
are also observed by the reflectance spectrum up to 13 
eV at room temperature for UO2 The energy-loss 
L(w) spectrum can demonstrate not only one-particle ex- 
citations but also collective excitations. The maxima at 
around 15.4 eV as showed in Fig. 2(c) indicates the plas- 
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FIG. 3: (a) The projected orbital-resolved partial DOS for 
Pu 6d, Pu 5/, and O 2p orbitals in antiferromagnetic PUO2. 
The Fermi level is set to zero, (b) The dynamical dielectric 
function e(u) — £i(<^») + ie2(u}) as a function of the photon 
energy ui for PUO2. The black and green lines represent our 
calculated real and imaginary parts of dielectric function e(w), 
respectively. 



mon resonance, which is qualitatively consistent with the 
experimental value of 14 eV l 20] . As showed in Fig. 1 (b) , 
at about 11.7 eV the real part e\ becomes zero, arriving 
at the minima around 12.1 eV and then approaches zero 
at about 14 eV. As Schoenes pointed out, the energy at 
which £1 (uj) crosses the zero line with a positive slop gives 
the plasmon excitation energy (2(j| . 



B. electronic structure and optical properties of 
PuO a 



Due to Pu unique position of its 5/ electrons between 
localized and dclocalizcd states among actinide series, Pu 
metal and plutonium-based oxides have more complex 
properties than other actinides. For example, metallic 
Pu has six different phase under different temperatures 
and pressures [23| . Pu02 as an important actinide diox- 
ide has extensive applications in nuclear reactor fuel and 
long-term storage of surplus plutonium. Therefore, the 
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study of optical properties for PuC>2 is also necessary and 
interesting. However, no experimental results of optical 
properties for PuC>2 in literature are available. Recently, 
Butterfield et al. studied the photoemission behavior of 
surface oxides of <5-plutonium and they observed that two 
peaks characterized by Pu 5/ and O 2p orbitals are dom- 
inant in PuC>2 and PU2O3 [2J|. For PuC>2, the two peaks 
observed are located at approximately 2.5 and 4.6 eV 
[24j . and our calculated DOS showed in Fig. 3(a) also 
present two similar peaks, i.e., a strong peak at about 
1.6 eV and a weaker one at 3.7 eV. Overall, these fea- 
tures are well reflected in our PDOS showed in Fig. 3(a) 
compared to experimental observations. As for the unoc- 
cupied 6d states, no experimental data can be obtained. 
Our calculated unoccupied 6d states begin at about 5 
eV. Considering the O 2p peak at —3.7 eV, we suggest 
the p — > d transitions occur at larger than 9 eV. 

Our calculated macroscopic dielectric constant £00 for 
Pu02 are also collected in Table I. The present e m is 6.21, 
whereas, no experimental value is unavailable at present. 
Our calculated imaginary part £2^) and real part £\{uj) 
of the complex dielectric function e(u>) are showed in Fig. 
3(b). For £2(0;), four main peaks locate at 2.8, 5.1, 7.5, 
and 10.5 eV. According to our PDOS calculation showed 
in Fig. 3(a), we attribute the peak below 6 eV to be intra 
5/ transitions, and the last two to be / — > d and p — > d 
transitions, respectively. The two similar peaks at 7 and 
10 eV are also obtained by Jomard et al. using ab initio 
calculations [25l |. 

Other related optical spectra for Pu02 are showed in 
Fig. 4. For reflectivity R(uj) spectrum, there are four 
peaks at 2.8, 5.0, 7.5, and 10.6 eV. Similarly, four peaks 
at 3.0, 5.5, 7.6, 10.6 eV are also observed in the adsorp- 
tion coefficient I(u>) spectrum. The origin of these peaks 
can also be explained according to the structure displayed 
in the imaginary part £2(w) of the dielectric function. It 
is evident that the plasmon excitation occurs at 16.0 eV, 
which is similar to the case of UO2 at 15.4 eV as men- 
tioned above. 



IV. SUMMARY 

In summary, we have performed a detailed investiga- 
tion of the electronic structure and optical spectra of ac- 
tinide dioxides UO2 and Pu02 using first-principle meth- 
ods. For UO2, our calculated projected orbital-resolved 
PDOS for U 5/ and O 2p orbitals in the valence region 
agree well with the experimental photoemission observa- 
tion. As for the unoccupied states, our calculated p-d 



gap is 5.8 eV, similar to the experimental BIS value of 
5.0±0.4 eV. The calculated insulating band gap E g and 
macroscopic static dielectric constants £oo for UO2 are 
2.3 eV and 5.53, respectively, which are also in good agre- 
ment with the experimental values of about 2.1 eV and 
5.1. The main features in spectra for UO2 are also well 
reproduced by our calculated dynamical dielectric func- 
tion e{oj) compared to the experimental observation. For 
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FIG. 4: Calculated optical spectra for PuCh, (a) the reflec- 
tivity R(oj), (b) adsorption coefficient I(oS), (c) energy-loss 
L(ui), and (d) refractive index n(uo). 

PUO2, the two main peaks characterized by Pu 5/ and O 
2p orbitals in valence bands are evidenced in our calcu- 
lated PDOS, which accords well with the photoemission 
results. The calculated macroscopic static dielectric con- 
stants £oo is 6.21. The related optical spectra for Pu02 
are also obtained by calculating the dynamical dielectric 
function. The / — > d and p — > d transitions are found 
to occurr at 7.5 and 10.5 eV, respectively. Considering 
the satisfactory optical description for UO2 compared to 
experiments, we expect that these results for Pu02 are 
also reasonable and therefore can provide a useful refer- 
ence for future experimental measurement. 
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